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RHIC and LHC. 
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I. INTRODUCTION 

The theory of strong interactions, QCD, incorporates 
the non-perturbative phenomena of confinement and chi- 
ral symmetry breaking. These phenomena are expected 
to disappear at high temperatures and/or at large chem- 
ical potential, where QCD can be analysed reliably be- 
cause its effective coupling becomes weak. How the phe- 
nomena arise in low temperature and low chemical poten- 
tial region is crucial to our understanding of the QCD dy- 
namics, ranging from hadronization in the early universe 
to the interior of neutron stars. A lot of effort has been 
devoted towards this goal, experimentally through the 
study of heavy ion collisions, theoretically through con- 
struction of phcnomcnological models (since we have not 
been able to solve QCD accurately), and computation- 
ally through simulations of lattice QCD. The results are 
often summarized as the phase structure of QCD in the 
temperature-chemical potential space. Recent reviews of 
the subject are available in Rcfs. . 

Consider QCD with N colors and Nf degenerate quark 
flavors of mass m, at temperature T and quark chemi- 
cal potential fj,. The phase structure of this theory is 
depicted in Figfj] in a schematic manner. (More details 
are needed to describe the real world QCD. Complexi- 
ties arising from unequal quark masses and color super- 
conductor phases at large chemical potential are omitted 
here.) Certain boundaries of the phase structure are bet- 
ter understood than the interior region, because of the 
exact symmetries present there. Explicitly: 

(1) For m — oo, the pure gauge theory has a finite tem- 
perature deconfinement phase transition, governed by 
the breaking of the global Zjy center symmetry of the 
Polyakov loop. This transition is of first order for N > 3. 

(2) For m = = jU, the theory has a finite temperature 
chiral phase transition, governed by the restoration of 
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FIG. 1: Schematic description of the phase structure of QCD 
in the m — T — fj, space. First-order transition surfaces are 
shown shaded, and critical lines are shown dotted. Value 
of m corresponding to the real world QCD is indicated by 
an arrow. Color superconductor phases occurring at large 
chemical potential are omitted. 



the flavor SU(N f ) v symmetry to SU(N f ) L ® SU(N f ) R . 
This transition is of first order for Nf > 3. 
(3) For m = = T, the theory has a baryon condensation 
phase transition, where the vacuum structure changes 
from ("i/>V) 7^ to (ip'ijj) 7^ 0. This transition occurs 
roughly when \x equals the constituent quark mass, and 
is also of first order. 

First-order phase transitions are generically stable 
against small changes of symmetry breaking perturba- 
tions. So the above three phase transitions extend in- 
wards, to varying extent, from the boundaries of the 
phase structure. The behavior expected from phe- 
nomenological models and numerical calculations is that 
the first-order transition surfaces end in critical lines, as 
shown in FigfT] Specifically, for the physical values of 
the quark masses, there is no phase transition as T is 
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varied, unless /i is sufficiently large. In the analytic cross- 
over region connecting high and low temperature regions, 
the three nearby transitions mentioned above cause vari- 
ous QCD properties to change rapidly. In general, which 
property is affected how much by which transition de- 
pends on the property concerned. To verify all the the- 
oretical expectations, therefore, it is important to con- 
struct appropriate experimental observables that high- 
light the individual aspects of deconfinement, chiral sym- 
metry restoration and baryon condensation, when QCD 
is tested at high temperature and/or at large chemical 
potential. 

In this article, I focus on the signals in heavy ion col- 
lision experiments related to the deconfinement phase 
transition. The collisions produce a high energy non- 
equilibrium state of QCD matter, which subsequently ex- 
pands and cools. The experimental signatures observed 
at RHIC and LHC demonstrate that the initial fireball 
first (quasi-) equilibrates to a quark-gluon plasma (QGP), 
which then hadronizes and decays. The detected quan- 
tities are multiplicities and distributions of various types 
of hadrons. Since a sizeable fraction of the colliding nu- 
clei continues unscattered along the beam direction, the 
signals coming from the fireball are cleanly observable 
only in sufficiently transverse directions, say for pseudo- 
rapidity |ry| < r? m ax- These signals can be broadly sep- 
arated into two types: (a) those describing the initial 
non-equilibrium fireball state and its approach to equi- 
librium, and (b) those characterizing the properties of 
an equilibrated QGP. The former consist of particles that 
can manage to pass through the QGP, e.g. direct photons 
and leptons, heavy quark jets, high px jets and the ellip- 
tic flow. The latter consist of distributions of moderate 
Pt hadrons produced close to the surface of the fireball, 
and they are the ones I want to relate to properties of 
the QGP. 

The QGP is a strongly interacting medium, while the 
hadrons resulting from its decay do not interact much 
with each other. The hadronization process thus goes 
through a stage where the scattering events among the 
components emerging from the fireball drop rapidly. Af- 
terwards, the hadrons essentially propagate radially out- 
ward without scattering, although the unstable ones de- 
cay. From the observed detector signals, it is possible 
to backtrack the distributions of hadrons to this "kinetic 
freeze-out" stage of the fireball. This stage is analogous 
to the "last scattering surface" in the evolution of the cos- 
mic microwave background radiation (CMBR), and sev- 
eral methodologies developed to study the fluctuations 
there are useful in its analysis. 

It should be kept in mind that the QGP is strongly 
coupled in the cross-over region, and the decoupling 
of hadrons goes through several stages instead of be- 
ing instantaneous Hadronization of the (quasi- 
equilibrated QGP takes place at energy density e cr ~ 
lGeV/fm 3 and T cr ~ 175MeV. Shortly thereafter, in- 
elastic scattering of hadrons stops, resulting in "chemical 
freeze-out" when T c h G m — 170MeV. Upon further expan- 



sion, elastic and resonant scattering (mediated largely 
by pions) ceases, producing "kinetic freeze-out" when 
Tki n ~ 120MeV. The scatterings contribute significantly 
to the thermalization of the hadron momenta, but as- 
suming that the hadronic medium has low diffusion, we 
can see through them to the patterns in the QGP. 

The distribution of hadrons at the kinetic freeze-out 
stage has often been modeled as a thermalized hadron 
resonance gas. Such thermal models fit the observed 
multiplicities of various hadrons, but miss all the multi- 
particle correlations between hadrons. The strongly in- 
teracting QGP should have many correlations, and more 
accurate models are therefore needed to relate them to 
observable signals. In what follows, I describe how a flux 
tube model providing a physical picture of the deconfine- 
ment process predicts specific two-particle correlations 
in the distribution of hadrons, going beyond the single 
particle multiplicities. I first present a summary of the 
model, and then discuss the baryon number correlations 
predicted by it. 

II. THE FLUX TUBE MODEL [f| 
A. Phenomenology 

The flux tube model of QCD is motivated by the dual 
superconductor description of linear color confinement 
[9|, where condensation of color magnetic charges re- 
stricts color-electric fields to vortex-like configurations. 
Although an exact derivation of the model from QCD 
has not been found, the model describes strong coupling 
expansions in lattice QCD and has been phenomenolog- 
ically quite successful. A characteristic property of the 
flux tube is its energy per unit length, i.e. the string ten- 
sion a. Other than that, the flux tube has a finite width 
w and a persistence length a (arising from the tube's 
stiffness so that the flux tube has to go a certain dis- 
tance before it can freely reorient itself), both assumed 
to be of order k^ CD . 

The flux tubes have to obey the constraint of Gauss's 
law. So they terminate only on quarks, and interact only 
at TV-point vertices. These two features represent the 
invariant tensors 5 a b and e a f, c (in case of SU(3)) used 
to describe the meson and the baryon wavefunctions. 
Other multi-quark hadron states are phenomenologically 
not prominent, except for multi-nucleon nuclei, and so all 
other interactions among the flux tubes are ignored in the 
model Note that a glueball would be represented by 
a closed flux tube loop in this description. 

The finite temperature behavior of the model is gov- 
erned by the competition between the energy and the 
entropy of the flux tube configurations. At low tempera- 
tures, energy wins, keeping down the total length of the 
flux tubes. At high temperatures, entropy dominates, 
producing elaborate structures of the flux tubes all over 
the space. 

First consider the situation for the pure gauge the- 
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FIG. 2: Possible flux tube configurations connecting a static 
quark-antiquark pair, as the temperature is increased (from 
top to bottom), and when baryonic vertices are included (from 
left to right). 



ory with only static color sources, i.e. m — oo. As the 
temperature is increased the flux tubes oscillate more, 
and also produce more vertices. Some of the possible 
configurations are shown in Fig(2] In absence of ver- 
tices, corresponding to the SU(2) gauge theory, there is 
a second-order deconfinement phase transition when the 
flux tube length diverges and the quark-antiquark pair 
loses information about each-other's position. In pres- 
ence of vertices, corresponding to SU(N > 2) gauge the- 
ories, the flux tubes can percolate the space in a network 
before their length can diverge. That also allows the 
quark-antiquark pair (hooked on to the network) to lose 
information about each-other's position, but produces a 
first-order deconfinement phase transition. 

When finite mass quarks are included in the model, 
they can break the flux tubes by quark-antiquark pair 
production from the vacuum. Baryon number produc- 
tion at finite chemical potential can also break the flux 
tubes. These possibilities are illustrated in Fig|3J With 
the flux tube network breaking up, the strength of the de- 
confinement phase transition weakens as the quark mass 
is lowered from m — oo and/or as the chemical potential 
is increased. Numerical estimates show that for Nf — 3 
QCD at small chemical potentials, the first-order decon- 
finement phase transition ends in a critical line around 
m = 1.5GeV, as depicted in FigJTJ Although the cross- 
over region does not have any sharp behavior of the 
deconfinement phase transition, we can still investigate 
whether the percolating flux tube scenario suggests any 
detectable signal there. 



B. Quantitative Formulation 

It is straightforward to formulate the model on a lattice 
with spacing a. The flux tubes live on the links, while 
the quarks and vertices live on the sites. The variables 
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FIG. 3: A color-electric flux tube can break when dynamical 
quarks are included in the theory, (a) A flux tube produced 
by static color sources, (b) Its breaking by a quark-antiquark 
pair appearing from the vacuum, (c) Its breaking by a baryon 
appearing from the vacuum at finite chemical potential. 



take values 0, ±1, depending on the direction of the flux 
(for N > 2), as shown in FigJU The total energy of a 
flux tube configuration is 



E = aaJ2\ 



m J2\pu\ + u £m > (!) 



where v denotes the energy cost of an iV-point vertex, 
and / sums over the 2Nf spin and flavor quark degrees 
of freedom. For a cubic lattice, the constraint of Gauss's 
law at every site becomes 



)-^2p^f + N ^ = a ^ = °- ( 2 ) 



When Gauss's law is applied globally, the contribution of 
the flux variables cancels. We then have two ways of 
determining the total baryon number B, either from the 
quark variables pij or from the vertex variables qf. 
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The grand canonical partition function for the system 
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FIG. 4: The link and site variables for the flux tube model. 
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with Eq.Q implying that quark chemical potential \x is 
equivalent to vertex chemical potential Nfi. The grand 
canonical partition function is fully factorized by express- 
ing the Kronecker delta constraint at every site as 
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The sum over the variables rii^, Pij, qi can then be 
carried out explicitly, resulting in 
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The equivalence of quark and vertex chemical potentials, 
following from Eq.Q, makes Z[T, fi] invariant under the 
global symmetry transformation, 9i — > 9i + e, for any 
complex value of e. 

By extending the allowed values of the variables 



Pi j, qi to all integers llj], the grand canonical par- 



tition function can be converted to [1J 
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The systems represented by Eq.© and Eq.© are both 
in the universality class of the XY spin model in the pres- 
ence of an ordinary magnetic field as well as a Z(N) sym- 
metric magnetic field. Their phase structures are similar, 
and their couplings are related by 
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In quantitative analysis, Eq.([7|) has the convenience that 
its algebraic structure is suitable for treatment using well- 
established methods in statistical mechanics. 

The physical meaning of the site variables 9i can be 
uncovered by looking at the free energy of a static color 
charge in the system. Introduction of a quark at site 
j modifies the Gauss's law constraint there as 5 a h o — > 
5 a .-i. That makes the free energy of a static quark 



exp(-F g /T) = (expH^)) 



(9) 



In the finite temperature gauge field theory language, 
this free energy is given by the expectation value of the 
Polyakov loop Pj at the site j. So we arrive at the cor- 
respondence that 6i represents the phase of the Polyakov 
loop Pi, and the flux tube description of deconfinement 
is dual to the familiar Polyakov loop description of de- 
confinement [l3| . The important advantage of the flux 



tube description is that it provides a position space vi- 
sual representation of what happens as the temperature 
is varied in QCD. 

Another advantage of the flux tube formulation is that 
it can be numerically simulated without any fermion sign 
problem at finite chemical potential [14l - fl7j , since Eq. fU) 
involves only real positive weights while Eq.([7]) necessi- 
tates complex weights. 



C. Limitations 

The flux tube model described above does not incorpo- 
rate the chiral nature of the quark degrees of freedom. So 
it cannot properly address the features in the QGP cross- 
over region that can be considered consequences of the 
nearby chiral phase transition. Attempts have been made 
to combine Polyakov loop and Nambu- Jona-Lasinio mod- 
els [l8l - l20| . to study chiral properties of the QGP. They 
have treated the Polyakov loop only as a mean field with- 
out its full dynamics, but have still provided reasonable 
fits to static QGP properties extracted from lattice QCD 
simulations. They need to be extended by including the 
kinetic term for the Polyakov loop dynamics, e.g. J =/= 
in Eq.flT]), to obtain a good estimate of the fluctuations. 
One of the predictions of chiral symmetry restoration at 
high temperature is an increased production of baryons 
relative to mesons, due to a decrease in the constituent 
quark mass. The same enhancement is also predicted 
by the increase in the number of flux tube vertices at 
high temperature. This coincidence makes the predic- 
tion robust, and it has indeed been observed (2lj . but an 
estimate of its magnitude requires accurate treatment of 
both confinement and chiral dynamics. That is an exer- 
cise for the future. 

Without the chiral properties of quarks, the flux tube 
model is also unable to say much about the baryon 
condensation phase transition. On the other hand, in- 
fluence of a finite chemical potential on the deconfine- 
ment phase transition can be inferred from the be- 
havior of the grand canonical partition function. The 
particle-O-antiparticle symmetry of QCD corresponds to 
invariance of Z[T, p] under simultaneous sign flips of the 
variables n,i^, Pij, qi- The partition function is there- 
fore an even function of /i, i.e. Z[T, fi] = Z[T, —fj]. Con- 
vexity of the exponential function then implies 

Z[T,p]>Z[T,n = 0], F[T,n]<F[T,/i = 0]. (10) 

F\T, fi] is continuous across the coexistence surface, and 
dF = —SdT — Ad/i at constant volume. So for positive 
latent heat, it follows that introduction of a chemical 
potential decreases the transition temperature. 

Furthermore, as described earlier, a non-zero chemi- 
cal potential weakens the signal from the deconfinement 
phase transition. On the other hand, it enhances the sig- 
nal from the baryon condensation phase transition, ob- 
servable in baryon number susceptibilities. Baryon con- 
densation needs to be treated as a site percolation prob- 
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lem and not as bond percolation by flux tubes. That will 
not be addressed here. 



III. BARYON NUMBER CORRELATION 
SIGNALS 

Let us consider the flux tube scenario of what happens 
in the heavy ion collision experiments as the fireball of 
the QGP expands and cools. Even though the cross- 
over region relevant to experiments does not possess a 
single percolating flux tube network, it can still contain 
many finite clusters of flux tubes. Because the flux is 
directed, an obvious feature of every such cluster is that 
any neighbor of a vertex is an anti- vertex and vice versa 
(see Fig. 2(d)). During the evolution of the fireball, the 
total baryon number is conserved :22j , and so the vertices 
can only be locally pair-produced or pair-annihilated. As 
the QGP hadronizes, the flux tube clusters start break- 
ing up. After the chemical freeze-out stage, there is no 
more production or annihilation of vertices; every ver- 
tex ends up in a baryon and every anti- vertex ends up 
in an antibaryon. In the absence of subsequent large- 
scale diffusion, the radial propagation of (anti)baryons 
preserves the geometric pattern of (anti)vertices present 
at the chemical freeze-out stage. As a result, the an- 
gular positions of the (anti)baryons seen in the detec- 
tor can be backtracked to the angular positions of the 
(anti)vertices on the surface of the fireball at the chemi- 
cal freeze-out stage. This pattern of vertices on the sur- 
face of the fireball can then be analysed for correlations 
and fluctuations, using techniques similar to those used 
to analyse the temperature fluctuations in the cosmic mi- 
crowave background radiation (23j . 

The above description is a simple consequence of pick- 
ing the right variables to visualize confinement dynam- 
ics in QCD. Although quantitative estimates of the 
(anti)vertex distributions would require numerical sim- 
ulations, there are certain qualitative features that can 
be gleaned with much less effort. The fact that heavy ion 
collisions produce a sizeable number of antibaryons [2lj |. 
from an initial state that has none, means that (a) a good 
number of baryonic and antibaryonic vertices are pro- 
duced in the fireball, and (b) fragmentation of flux tubes 
during hadronization is more likely than their shrinking 
and coalescing that would annihilate vertices with anti- 
vertices. Our aim is to look for a specific pattern in the 
distribution of the produced (anti)vertices. Note that the 
flux tube scenario is unable to say anything regarding the 
energy and the momenta of the hadrons that emerge from 
the fireball. 



A. Experimental Data Parametrization 

There are several differences in the type of data gath- 
ered from CMBR and from heavy ion collision experi- 
ments, which need to be kept in mind when studying the 



two with similar techniques. 

(1) The CMBR is a single high statistics event where ho- 
mogeneity and isotropy of the universe allow accurate de- 
termination of the distribution parameters, while heavy 
ion collisions are multiple modest statistics events where 
ensemble averaging improves the accuracy of the distri- 
bution parameters. 

(2) The temperature data of CMBR are real numbers, 
while the baryon number data of heavy ion collisions are 
integers that need to be binned and smeared. 

(3) The average temperature and size of the "last scat- 
tering surface" are well-determined for CMBR, while the 
number of participating nucleons and the fireball volume 
fluctuate considerably in case of heavy ion collisions and 
are not accurately determined. 

(4) The CMBR data cover the full 4n solid angle and 
hence can be easily parametrized using the orthogo- 
nal basis of the spherical harmonics, while the heavy 
ion collision data arc restricted to sufficiently transverse 
directions only and a suitable orthogonal basis for its 
parametrization has to be found. 

Let us assume that the available data cover 6 m < 9 < 
Tt—9 m , and label the angular distribution by the unit vec- 
tor n. The heavy ion collision experiments have an axial 
symmetry around the beam axis as well as the reflec- 
tion (or parity) symmetry 9 -H- tt — 9. These symmetries 
are sufficient to prevent the average baryon number flow 
in any direction, and can be used to improve the sta- 
tistical accuracy of the correlations in the data. Fourier 
expansion helps in search for correlations by orthogonal 
separation of scales, and the baryon number distribution 
can be partially expanded as 

&(n) = &(M) = 4=E E W e "" - (n) 

V Z7T — 1 

v a— ± ra— — oo 

Here b± (6) = ±6±(tt-0), reality of b(n) implies V£{6) = 
b a L m (9), and (b(h)) = (9)/V2tt. In a sense, the positive 
and the negative parity terms correspond to all the even 
and all the odd values of respectively, of the spherical 
harmonics expansion. Since the hadrons emerging from 
the fireball in longitudinal directions go undetected, the 
total and detected baryon numbers of a collision event 
are (note that tanh?7 = cos 9): 

B tot = V2^ J d{cos9) b+(9) , (12) 

/cos 9 m 
d(cos9) b+(9) . (13) 
- cos 9 m 

The two-point baryon number correlation function is 
the ensemble average (b(h)b{n')) . Because of the axial 
symmetry, it depends only on the difference of the az- 
imuth angles <f> — </>', and because of the reflection symme- 
try the product of the parities aa' has to be one. There- 
fore, 

{b° m {9)b°: m ,{9')) = (b° m (9)b° m r(9>)) = 5 aa ,5 mm ,C° m {0,e') , 

(14) 
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and the two-point correlation function becomes 

1 oo 

(b(h)b(n')) = -J2 E C£(M') e*"^-*') . (15) 

(T=± 771 — — OO 

Inverting this relation, we obtain 

C£(M') = ^fj^fj^' e m ^'-V(b°(h)b°(n')) , 

&±(n) = i(6(e,0)±6(7r-e,^)) . (16) 

The coefficient functions C^(6,9') are real and symmet- 
ric. They contain all the information about the two- 
point correlation function. The non-trivial correlations 
are given by the connected contributions, 

(b(h)b(n'))c = (b(h)b(n')) ~ (b(n)) (b(ti)) , (17) 

[C^(e,6% = C? n (9,9') - 5 m0 S a+ b+(0) b+{6') . (18) 

In practice, depending upon the resolution available in 
the data, (a) the expansion can be truncated at a suitable 
value of |m|, and (b) the ^-dependence can be further 
subdivided into smaller bins. 

Similar parametrizations of the data can be carried 
out for distributions of other conserved quantities also, 
e.g. the electric charge or the strangeness. As a matter 
of fact, experimentally observed correlations for different 
conserved quantities can be compared with each other, 
to illustrate different aspects of the QGP dynamics. The 
flux tube model predictions are the cleanest in case of the 
baryon number, and illustrate the deconfinement mecha- 
nism. Different models with appropriate features would 
be required to understand the dynamics of other con- 
served quantities. 

It should be noted that the theoretical estimates of the 
correlations are obtained assuming that the volume of the 
QGP is sufficiently large. On the other hand, the number 
of participating nucleons in the collisions, and hence the 
volume of the fireball, depends on the centrality of the 
heavy ion collisions. So to determine as well as to con- 
trol the consequent systematic effects, it is necessary to 
determine the correlation functions in different centrality 
ranges separately. 



B. Theoretical Expectations 

Theoretically, the equilibrium correlations between 
flux tube vertices are straightforward to calculate in the 
three-dimensional position space, which can then be pro- 
jected onto the observed surface of the fireball. An al- 
ternating neighbor pattern of vertices and anti-vertices, 
expected in a percolating flux tube network, would give 
rise to two-point baryon number correlations similar to 
the two-point charge correlations in ionic liquids (with 



ions of comparable size) such as Cs + Cl~. The key com- 
mon ingredient in the two cases is the hard-core repulsion 
between the objects involved. 

For discrete objects located at r*; , a convenient descrip- 
tion of the position space correlations is in terms of the 
pair distribution function 

where p(r) = (^ S(r — fi)) is the local density of the 
objects. In homogeneous and isotropic fluids, p is inde- 
pendent of the position and g depends only on \r — f^\. 
Choosing r* = 0, then we have 

P9(r) = (j2 6 (r-n)) , (20) 

Interactions fade away at long distances, and so asymp- 
totically g(r oo) = 1. For objects with no correlations, 
e.g. an ideal gas, g(r) — 1. For objects with hard-core re- 
pulsion, g(0) = 0, and beyond the hard core g(r) tends to 
its asymptotic value exhibiting damped oscillations [24j . 
as illustrated in FigJS] In particular, the distance scale 
of the oscillations is determined by the inter-object sep- 
aration (the first peak corresponds to the likely nearest- 
neighbor separation, the second peak corresponds to the 
likely next-nearest-neighbor separation, and so on), and 
the amplitude of the oscillations is determined by how 
tightly the objects are packed together. 

To detect the correlations between vertices and anti- 
vertices, it is useful to compare the pair distribution func- 
tions for \q(r)\ |g , (f a ')| and q(r) q{f), g\ v \ and g v respec- 
tively, on the same data sets. Comparison of g\ v i and 
g v keeps under control the effects arising from variation 
in the number of participating nucleons in the heavy ion 
collisions. Both follow the form of Eq. fp~9|) . with the sum 
running over vertices and anti- vertices. The former omits 
any vertex value signs and uses the density pi v i, while the 
latter includes the vertex values, i.e. qiqj on the right- 
hand side of Eq. (ll9p and qi in obtaining the density p v 
(note that p v < p\ v \). 

For a percolating flux tube network, the function 
g\ v \(r) should behave similar to that for hard-core ob- 
jects, with the nearest-neighbor separation of the order of 
the baryon size. Similar behavior would be expected from 
several other fluid models also, so g\ v \(r) can be treated 
as a model-insensitive reference function. In contrast, 
in a percolating flux tube network, successive neighbors 
contribute with opposite signs to g v (r). The resultant 
pair distribution function would then have oscillations of 
the form sketched in Fig[51 For a system with little or 
no correlations between vertices, e.g. single particle ther- 
mal distribution models, the probability of occurrence of 
a vertex or an anti- vertex at any location is proportional 
to its overall density. Contributions to the right-hand 
side of Eq. (fTTH) would then have the same density factors 
as on the left-hand side, and g v (r) would behave the same 
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FIG. 5: Schematic representation of the pair distribution 
functions g\ v \(r) (thick line) and g v (r) (thin line). The former 
is similar to that for objects with hard-core repulsion. The 
latter is for a percolating flux tube network where vertices 
and anti- vertices alternate. 



way as g\ v \(r). Thus <?„(r) is sensitive to the correlations 
present between vertices and anti- vertices, and the con- 
trast between <7|„|(r) and <7„(r) can be used as a measure 
of these correlations. 

The percolating flux tube network in the high temper- 
ature QGP is maximally connected for m = oo, with N 
nearest neighbors for each vertex and anti-vertex. The 
network starts breaking up as the quark mass is low- 
ered, more so if the chemical potential is also increased. 
Even individual baryons that are not part of any flux 
tube cluster can appear. Such a break up reduces the 
correlations between neighbors, smoothening the oscilla- 
tions of the pair distribution function, analogous to what 
happens when a liquid is heated up. For the physical 
values of the quark masses, there is no deconfinement 
phase transition and no single percolating flux tube net- 
work is expected in the QGP. Still, if large enough flux 
tube clusters survive in the cross-over region, they would 
contribute to the contrast between g v (r) and g\ v \ (r). The 
strength of the contrast, therefore, measures the extent 
to which the nearby deconfinement phase transition influ- 
ences the properties of the QGP in the cross-over region. 

In general, the first peak of the pair distribution func- 
tion is the most informative about the fluid properties. 
Its location provides an estimate of the inter-object sepa- 
ration, and the area under it is a measure of the number 
of nearest neighbors. In the present case, it is also the 
place where the contrast between g v (r) and g\ v \(r) is the 
maximum, and so it is the best feature for identifying the 
correlations between vertices and anti- vertices. Note that 



for specific values of the parameters in Eq.©, numerical 
simulations of the flux tube model can obtain the func- 
tions g v (r) and g\ v \(r) for an equilibrated homogeneous 
and isotropic QGP. 



C. Connecting Theory to Experiment 

Projection of the three-dimensional pair distribution 
function onto the surface of the fireball smears the os- 
cillatory structure of g(r). The resultant angular distri- 
bution can be expressed in terms of a smearing function 
as 



S(a,r) g(r) dr , / S(a,r) da = 1 



w(a) - 

J r min J o 

.(21) 

Here the integral over r has to be sufficiently restricted 
so as not to lose the oscillatory signal. Since the corre- 
lations that we are interested in are of short range, it is 
convenient to remove the uniform part of the distribu- 
tion function, e.g. look at g(r) — 1 which describes the 
connected part. In the context of picking up the vertex 
correlation signal, it is even better to look at the differ- 
ence <7„(r) — g\ v \(r), which vanishes for both r — > and 
r — > oo. 

In a specific setting, the smearing can be performed 
with sufficient accuracy numerically. Approximate an- 
alytic result can be obtained, nevertheless, under the 
assumption that the inter-vertex separation r is much 
smaller than the radius of the fireball R. Let a be the 
position of the center of the pair with respect to the cen- 
ter of the fireball, and let /? be the angle between a and f. 
Then the pair separation r projects to angular separation 
a = (r/a) sin/3. Taking expectation value over j3 and a, 
for a homogeneous and isotropic distribution of pairs in 
the fireball, the smearing function becomes: 



S(a,r) = 



R 3a 2 da f n sin /3 dp 



Jo R 3 -hi 

f R 3a 2 da 



aa 



(a - - sin/3) (22) 



R 3 



97rr 3 /(16i? 3 a 4 ) :a>r/R 

(12A - 8 sin(2A) + sin(4A)) (23) 
: sin A = aR/r < 1 



:i2/V ;; 



This small angle approximation to S(a, r) should be rea- 
sonable for studying at least the first two peaks of g v (r), 
given that experimentally the radius of the fireball in cen- 
tral heavy ion collisions is ~ 6fm and the inter-nucleon 
separation in nuclear matter is ~ 2fm. 

When the rotational symmetry is a good approxima- 
tion, it is convenient to decompose the angular distribu- 
tion w(a) in terms of the orthogonal Legendre polyno- 
mials: 



w(a) = J2 C l 



1=0 



21 + 1 

47T 



Pi (cos a) 



8 



IV. SUMMARY AND OUTLOOK 



Ci = 2?r / d(cos a) w (a) Pi (cos a) . (24) 



An estimate of w(a) can then be compared to experimen- 
tal data using cos a = cos 9 cos 9' + sin 8 sin 9 1 cos(</>— </>'), 
but Ci cannot be extracted from the experimental data 
unless full angular coverage is available. With only ax- 
ial symmetry and limited angular coverage present, w(a) 
needs to be decomposed in terms of parity and azimuthal 
Fourier components as in Eqs. (|15ll6[) . Specifically, that 
avoids systematic errors in the construction of the con- 
nected distribution w c (a). The addition theorem for the 
associated Legendre polynomials, 

i 

P ; (cosa)= (-l) m J Pr(cos6»)Pf m (cos6»')e jm( ^*' ) , 

m——l 

(25) 

then relates the coefficient functions according to 
C'(M') = (-1)™ (26) 

OO 

x £' c M^) p n cos0 ) p r m (cos0'), 

where the sum is restricted to even(odd) values of I — |m| 
for a = +(— ). The functions 0^(0,6') can be extracted 
from theoretical models as well as from the experimental 
data, and they can be restricted to suitable ranges of \m\ 
and 9 depending on the resolution available. 

There still remain gaps between the theoretical for- 
malism described here and the actual experimental data. 
The most prominent one is the fact that the detectors 
record only the charged hadrons. Thus protons and anti- 
protons are observed but neutrons and anti-neutrons are 
not. The requisite baryon number correlations can be 
extracted from the data only if the observed subset of 
(anti)protons provides a faithful characterization of the 
total baryon number distribution — ideally the two should 
be proportional. Moreover, corrections need to be es- 
timated due to only approximate equilibration of the 
fireball, non-uniformity of the QGP caused by the el- 
liptic flow, and baryon number diffusion subsequent to 
hadronization, all of which are likely to weaken the os- 
cillatory correlation signal. On the other hand, develop- 
ment of the hard baryon core during hadronization would 
enhance the oscillatory signal, compared to the softer 
QGP state at higher temperature. Despite these gaps, 
it is worthwhile to look for the two-point baryon num- 
ber correlations in the experimental data, as a charac- 
teristic signature of the deconfinement phase transition. 
Specifically, experimental determination of the contrast 
between g\ v \ and g v , and their departure from g(r) = 1 
corresponding to no correlations, would provide a model- 
independent characterization of the baryon- antibaryon 
correlations. 



The Polyakov loop is a widely used order parameter 
for understanding the finite temperature deconfinement 
phase transition in a gauge theory. Its dual description 
in terms of the flux tube model has the advantage that 
the gauge theory dynamics can be visualized in position 
space. This picture helps in connecting experimentally 
observable baryon number correlations with properties 
of the deconfinement phase transition. There is no fun- 
damental interaction associated with the baryon number. 
So the baryon number correlations have to arise from its 
precursors in the QGP, i.e. the pattern of the flux tube 
vertices expected from QCD. 

The experimental effort in studying heavy ion collisions 
has so far largely focused on single particle distributions. 
In this article, I have gone beyond that towards analysing 
multi-particle distributions. I have described how to ex- 
tract the two-point baryon number correlations from the 
experimental data, and have also predicted that it would 
contain an oscillatory signal, based on theoretically ex- 
pected QGP properties at the chemical freeze-out stage. 
The experimental data for moderate px baryons in cen- 
tral heavy ion collisions are the best suited for observa- 
tion of the predicted signal. The features to be quan- 
tified are the distance scale and the amplitude of the 
oscillatory correlations. The former is expected to be the 
inter-baryon separation (~ 2fm) , while the latter though 
uncertain at present would tell us a lot about how tightly 
or softly the QGP is packed. It is worth keeping in mind 
that the intense search for two-point correlations in the 
CMBR, theoretically expected but without any accurate 
prediction of its magnitude (inflation is a much more flex- 
ible theory than QCD), found the correlations at the level 
of AT ss I0 -5 T, and provided a major boost to our un- 
derstanding of cosmology. 

As a different application, the flux tube picture may 
also help in improving the baryon production estimates 
of hadronization models, by replacing the creation of di- 
quark pairs with the creation of baryonic vertex pairs. 
It can also be noted that the flux tube dynamics influ- 
ences diffusion in the QGP. Specifically, formation of a 
flux tube network would suppress diffusion. Both dif- 
fusion and viscosity characterize transport properties of 
a fluid, incorporating dissipative dynamics and entropy 
production. The former is convenient to use in position 
space, the latter in momentum space, and the two can 
be related in a kinetic theory framework. Low diffusion 
is not incompatible with the small viscosity to entropy 
density ratio in the high entropy density QGP. Quanti- 
tative analysis of diffusion in the QGP requires a hydro- 
dynamic framework, however, which is beyond the scope 
of the formalism discussed here. 
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